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Abstract 

Using a sequence of rational approximants and the large- A'c limit of QCD, 
we estimate the value of the low-energy constant Cgy which appears in the La- 
grangian of Chiral Perturbation Theory at 0{p^). 



The Chiral Lagrangian[Tl [2] organizes the physics of the strong interactions at low 
energy as an expansion in powers of momentum and masses of the hghtest pseudoscalar 
fields, which are the only ones explicitly present in this Lagrangian. Since all the 
heavier states of QCD are integrated out, their physics is encoded in a set of low-energy 
constants (LECs). These LECs are indispensable to make definite predictions in Chiral 
Perturbation Theory. There is already a relatively good knowledge of the value of most 
of the LECs which appear at O(p^) in the chiral expansion[2], [3]. At 0{p^), however, 
most of the (9(100) LECs are completely unknown. This note is concerned with the 
estimate of one of them, the LEC of 0{p^) CgrflO] appearing in the < VV — AA > 
two-point correlator in the chiral limit. 

The general strategy will consist in constructing a rational approximant to the rele- 
vant Green's function, i.e. < VV — AA >, from the coefficients of the chiral expansion 
and any other known properties of the full function. Once the rational approximant is 
known, upon reexpansion around = 0, higher order unknown coefficients of the chi- 
ral expansion may be predicted. If the rational approximant is a better description of 
the original function than the partial sums of the chiral expansion, one may expect this 
prediction to be reliable. For a brief review, where further references to the literature 
may be found, we refer to Weniger|[l]. 

Let us, therefore, consider the two-point functions of vector and axial-vector cur- 
rents in the chiral limit of QCD 

Kl^'i^) = ' I d'xe'^^J^^\x)Jl ^'^(0)) = {q,q. - g,.q') UyAq') , (1) 

with Jy{x) = d{x)'~f'^u{x) and = d{x)'j'^'j^u{x). As it is known, the difference 

Uv-Aiq'^) satisfies the unsubtracted dispersion relation given bj|l] 

1 df 1 

Hy.Aiq') ^ -{Jlv{q') - n^(g^)) = lim / ^ - Imnv._A(t) • (2) 

Z A^oo Jq t — q^ — le tx 

Since all LECs are defined in the chiral limit, the restriction of the function liy-Aiq^) 
to this limit entails no loss of generality. Even then, the analytic structure of Ilv-Aiq^) 
is very complicated, with a multiparticle cut starting at t > 0. A further simplification 
occurs in the large- A^^^ limit of QCD [6] in which the previous cut becomes suppressed 
and only single particle intermediate states are allowed. The function Iiy_A{q^) con- 
tains then an infinite set of isolated poles [9], and becomes meromorphic. In Ref. [7] 
it was emphasized that any truncation of this meromorphic function to a finite set of 
poles may be interpreted as a rational approximation to the original function. There 
are a number of reasons why the large- A'c limit of QCD is a sensible limit to take, in 
particular for Green's functions built out of the vector and axial- vector currents [SIE]- 

Therefore, in the large-A^c limit, the function g^ny_^(g^) has the following repre- 
sent at ioiJl 

oo 

q'^v-A{q') = fl + q'i:—^ (3) 

^The upper cutoff which is needed to render the dispersive integrais matheniaticaiiy weii defined 
can be sent to infinity provided it respects chiral symmetry [5] . 
^Multiplication by kills the pion pole at the origin. 
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where R labels resonance states and, assuming the existence of a mass gap, all masses 
are nonzero with the rho meson mass being the smallest one in the sum. In this way, 
defining = — g^, the expression is analytic at = 0, and allows the Taylor 
expansion |26j 




where /g is the pion decay constant (in the chiral limit), and Liq, Csj are the LECs at 
0{p^,p^) (respectively) in the corresponding Chiral Lagrangianj2| [TO]. On the other 
hand, the expansion of Ily.yi at infinity is 



-Q^Uv.A{-Q')^^7ia,'-^(l + 0{as\ogQ'))+... , (Q^ ^ oo) (5) 



where as is the QCD coupling constant and (V'V') is the quark condensate. Unlike 
the expansion around the origin, the existence of nonvanishing anomalous dimensions, 
even in the large- A'^c limit, gives rise to the logQ^ terms and, unlike (jlj), renders the 
expansion around infinity in ([5]) not analytic. 

The expansion (jlj) will be our starting point in this work. Although the LEG Liq 
is pretty well known [21 [TT]. this is not so for Cs7- It is therefore important to obtain a 
new determination of this LEG with its associated error. 

Given the meromorphic function Q'^IIv-a{—Q^) in the complex plane with an 
analytic expansion around the origin, as in (jlj), it is possible to construct a Pade 
Approximant (PA), P^^(Q^), as the ratio of two polynomials of degree M and A^ 
(respectively) such that its expansion in powers of matches that of the original 
function up to, and including, the term of (9(Q^*^^"'"^)). Since the function falls off 
at large as up to logarithms (see ([5])), we choose A^ = M + 2 in order to 
optimize the matching of the rational approximant at large to this behavioiE]. We 
emphasize, however, that this choice does not affect the properties of convergence of 
Pade Approximants, as described next. 

As M — »• oo, there is a theorem|T2] that ensures convergence of the sequence of 
PAs to the original meromorphic function, in any compact set in the complex 
plane except at a finite number of poles. Of course, where there is convergence, the 
PA may be considered an approximate resummation of the Taylor series around the 
origin. On the other hand, the set of points where there is no convergence certainly 
includes the position of the poles since not even the original function is defined there, 
but there may appear other artificial poles which have no counterpart in the original 
function. One would naively think that the presence of these artificial poles would 
cause a major distortion and completely spoil the rational approximation. However, 
one can actually show[12j that as the order of the Pade increases, i.e. as M grows, 
these artificial poles either move to infinity in the complex plane and decouple or they 
get "almost-canceled" by the appearance of nearby zeros. Although, in general, this 
cancelation is not complete, it is efficient enough to make the region of distortion of 
the artificial pole only of zero measure. This is why and how the Pade Approximation 

■^Due to the presence of logarithms in ([S]), however, this matching cannot be perfect. 




Q^Uv-Ai-Q'') ~ /o - 4 Lio - 8 Csr + . . . 



iQ' - 0) 



(4) 
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works. For an explicit example where all these properties come to play in the context 
of a Regge-inspired model, we refer to [12]. 

Theorem [12j is important because it teaches us useful information about the qual- 
itative behavior of how Fades approximate meromorphic functions. Regretfully, when 
asking more quantitative questions such as the rate of convergence, which is the first 
step towards an estimate of the error, such a theorem is only of limited practical im- 
portance. In practice, one can take a more useful approach towards an estimate of the 
error by studying the behavior of a set of successive rational approximants, as we will 
now explain. 

In order to be able to construct a sequence of rational approximants it is of course 
crucial to have enough number of inputs. Since FAs are constructed from the coeffi- 
cients of the Taylor expansion (j4j) one immediately faces an obvious difficulty. Since 
what one wishes is an estimate of Cgy, only the two coefficients /o and Liq may be used. 
With these two coefficients as input, the only FA vanishing at large is P°, but it 
falls off as which is too slow as compared to ([5]). Consequently, it is necessary to 
consider more general rational approximants than the standard FAs. 

In Ref. [TS] we saw that the 7r+ — tt^ mass difference in the chiral limit, which is 
given by 

o /"OO 

6M'^ ^ Ml, - Ml, = ^ / dQ' ny.A{-Q') , (6) 

47r /o Jo 

could be used as a further constraint in the construction of the FAs with very good 
numerical results. Together with /o and Liq one now has three inputs to construct the 
P2, which does match the power fall-off at large in ([5]). By simple re-expansion 
around = it is then possible to predict an estimate for the term of 0{Q'^) in (jl]). 
This prediction was checked against the exact value in the model in [ISj with very good 
results, and this encouraged us to do the same also for QCD. In the QCD case, using 
the valued 

/o = 0.086 ± 0.001 GeV , (Ref. [12]) 
SM^ = 4.5936 ± 0.0005 MeV , (Ref. [20]) (7) 
Lio(0.5 GeV)<Lio<Lio(l.l GeV) =^ Lw = (-5.13 ±0.6) x 10"^ , (Ref. [2l]) , 

we then obtained [18] an estimate of the 0{Q^) term in the expansion (jl]) which 
translates into the value Cgr = (5.4 ±1.6) x 10"'^ GeV~^. In the present work, we 
would like to reassess this value with a more complete analysis. 

There are different kinds of rational approximants closely related to the usual FAs 
which, although perhaps not so popular, are also very useful. Among those, we would 
like to stress the so-called Fade- Type Approximants (FTA)[22J, TfJ{Q'^), which are 
very useful when one has some knowledge of the spectrum of resonance masses in the 
original Green's function. They are defined as the ratio of two polynomials Qm{Q'^) 

^ Since nis decouples from Ft^ in the large- iVc limit, the value of /o is estimated in Eq. ([7|) by 
extracting the chiral corrections from F^r using SU{2) x SU (2) chiral perturbation theory, but doubling 
the error as compared to Ref. [TO] , 
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and T]\f{Q'^) (of degrees M and A^, respectively): 

where the polynomial in the denominator has its zeros preassigned precisely at the 
positions of the first resonance masses in the original Green's function ([3]), i.e. 

Tn{Q') = iQ' + M2)(Q2 + m2)...(Q2 + M^) , (9) 

and the polynomial Qm{Q'^) is defined so that the expansion of the PTA around = Q 
agrees with that of the original function up to terms of order M + 1, i.e. 

Tj^'(Q') ^h + hQ' + f2 Q' + ... + fM Q''' + 0(g'(*'+^)) . (10) 

Choosing N = M + 2 one optimizes the matching of "^^m+2{Q'^) expansion ^ 

at ^ cxD and, as in the case of PAs, this is a choice we will make. 

Both PTAs and PAs where studied in Ref. [IH] and the lessons which can be 
drawn from that model are the following. The model confirms that one may estimate 
the unknown LECs with these rational approximants where, in the case of PTAs, the 
physical masses were chosen in increasing order, i.e. Mi < M2 < M3... For instance, 
with the PTA '^M+2iQ'^) could see that one has a good prediction for the term of 
(9(Q^(*^+^)) in the low-Q^ expansion, which is the first one not used as input, with a 
precision which improves as the order of the approximant, M, increases. Furthermore, 
the accuracy obtained for the unknown coefficients of the Taylor expansion is very 
hierarchical: the accuracy obtained for the term 0{Q'^^^^~^^^) is better than that for 
the term of 0(Q2(Af+2))^ ^j^^^^ better than for the term 0(Q2{m+3))^ ^j^h a quick 
deterioration for higher-order terms. The case of PAs follows the same pattern. As to 
the description of the spectrum, we found that PAs also reproduced the values for the 
residues and masses in a hierarchical way: while the first masses and residues are well 
reproduced, the prediction quickly worsens so that the last pole and residue of the PA 
has no resemblance whatsoever with its physical counterpart. The same is true for the 
residues of a PTA (since the masses are fixed to be the physical ones by construction). 

Based on the above, one can envisage the following strategy for getting a sequence 
of estimates for the 0{p^) LEG Cs7- Assuming that the vector and axial-vector meson 
masses stay approximately the same in the large- A'c and chiral limits, one can use their 
values extracted from the PDG book^20j to construct several PTAs. We think that this 
assumption is reasonable for both limits. First, for the chiral limit, this is because the 
up and down quark masses are very small[T5]. Second, for the large- A'c limit, there is a 
non negligible amount of phenomenological evidence in favor of the rho meson being a 
qq state [16], [8]. Besides, the success in the spectroscopy of the quenched lattice results 
for the lightest vector mesons is also suggestive that l/N^ corrections may not be very 
large [ISjj. Therefore, we will use for the masses 

rup = 0.7759 ±0.0005 ,mp, = 1.459 ± 0.011, m^/. = 1.720 ± 0.020, m^/// = 1.880 ± 0.030 

= 1.230 ± 0.040 , m„'^ = 1.647 ± 0.022, (11) 



^Be that as it may, whether the assumption is correct or not will ultimately be judged by the final 
results obtained. 
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where all the numbers have been expressed in GeV. 

For instance, with only /q and the masses of the p and Oi, one can construct the 
PTA T0(g2) and predict the value for = (-4.32 ± 0.02) x lO'^, which is not bad 
when compared, e.g., with ([7]). The next term in the expansion gives the following value 
for Cs7 = (4.00 ± 0.09) x lO'^ GeV-^ which is similar to that obtained in [18] with 
the Pade • However, since this value for Cgy comes from the second unknown term 
in the expansion of T2(Q^) rather than the first, it is quoted here only for illustrative 
purposes and will not be included in our final estimate, in agreement with our previous 
discussion. Adding Lio and the p' mass to the previous set of inputs one can then 
construct TUQ"^), which produces Cgr = (5.13 ±0.26) x lO'^ GeV-^. The PTA can 
be constructed if one also uses the pion mass difference ([6]) and m^'^, yielding in this 
case Cs7 = (5.24 ±0.33) x 10"^ GeV"^. We find the stability of these predictions quite 
reassuring. 

A comment on the quoted error estimates is in order. These quoted errors are 
the result of the propagation of errors from the input via the montecarlo method [T7]. 
As such, they do not reflect the intrinsic systematic error due to the approximation 
itself which will be estimated, at the end, as the spread of values obtained with the 
sequence of different approximants. On the other hand, the propagation of the error 
from the input via the montecarlo method consists in the following. Taking each input 
in ([7j) and (fTTI) . we have constructed a sample of data with a gaussian probability 
distribution yielding as the average and standard deviation precisely the corresponding 
input value and its quoted error, respectively. For each member of this sample, the 
rational approximant is then constructed and, upon reexpansion, the LEG is obtained. 
The distribution of the different values for Cgy so obtained happens to be also gaussian 
to a very good approximation. Therefore it will have an average value X and a standard 
deviation Y which are then used to quote the result for Cgr as X ±Y. 

To be able to construct further rational approximants one needs an extra assump- 
tion. Although, as we have emphasized above, the residues of the heaviest poles in a 
rational approximant do not come out anywhere close to the corresponding physical 
decay constants, this is not true for the lightest ones. In particular, in Ref . [18] , it was 
seen that the value of the residue for the first pole in a PTA could reproduce the exact 
value in the model with very good precision if the order of the PTA was high enough 
and, more importantly, it was improving as the order of the PTA grows. Gonsequently, 
if we are willing to use the decay constant Fp, and perhaps also the F^j, one can go for 
the construction of higher PTAs. These two residues can be gotten from the decays 
p — > e~^e~ and ai tt'J, respectively, and their values arepB] 

= 0.156 ±0.001 =0.123 ±0.024 (12) 

in GeV units. 

For instance, using /q , Lio, SM^ and Fp, as well as the five masses rup, ma^.mpi .ma'^ 
and mp'/, one can construct the PTA T|. Upon expanding this approximant, one 
obtains the value Cgy = (5.78 ± 0.21) x 10"'^ GeV~^. Alternatively, one can also use 
/q,Lio, Fp and only the first four masses to construct a T| approximant, which is 
different from the other T4 considered above. The value obtained for Cst, i.e. C^-j = 
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inputs 




/o ; TUp, rua 


T3 


fo, Lio ; rrip, rria, nip' 


^2 (a) 


/o, Lio, 6M^ ; rup, rUa, m^', m^' 


^2(6) 


/o, Lio, Fp ; mp, ma, m^', m^> 


r|p3(a) 


fo, Lio, Fp, ; rup, nia, nip', m^' , m^" 


^3(6) 


fo, Lio, Fp, Fa ; nip, m^, m^, , m^" 


■"■6 


fo, Lio, Fp, Fa, SM^ ; nip, nia, nip', nia' , m^", m^'" 



Table 1: Set of inputs used for the construction of the different Fade Type Approxi- 
mants in the text. 

(6.00 ±0.15) X 10~^ GeV~^, is nevertheless very similar, which again brings confidence 
on the prediction. 

In this way we have constructed a variety of rational approximants which we have 
listed on Table 1, in increasing order of the degree in the denominator, together with 
the set of inputs used. We have gone all the way up until the Tg, with the six masses 
hsted on (fTTl). 

Figure 1 shows the prediction for the LEG Csj from the corresponding rational 
approximant shown on the abscissa, upon expansion around = 0. We also included 
our previous result obtained in Ref. [18] with the FA P2, but with the present mon- 
tecarlo method for the treatment of errors. As one can see, the stability of the result 
is quite striking. After averaging over all these points, we obtain as our final result in 
the large- A^^c limit, 

Cs7 = (5.7 ± 0.5) X 10-^ GeV-2. (13) 

The error in (fT3l) is mainly dominated by the error on the input for Lio in E<1- dZ]) and 
is rather insensitive to the errors on the other inputs. For instance, one could increase 
the error on fo to 5 MeV in Eq. ([7]), or the error on nip to 50 MeV in (fTTj) . or the error 
on 5Mn- to 0.5 MeV in ([7]), without falling out of the error band given in fll3p . 

For comparison, we also show in Fig. 1 the result of several previous estimates 
for this LEG. Reference [26] (shown as 'A') uses the residues in Eq. f|T2l) and the 
p and Ol physical masses to construct, in effect, what we could call the FTA ¥3 to 
Q^Ilv-A- The difference between this result and ours stems from the fact that this 
rational approximant falls off like a constant at large Q^, unlike Eq. Also, as we 
have already emphasized, the use of the physical decay constant Fa^ ffT^ in a rational 
approximant which has the oi as the heaviest pole is a potential source of error. 

Reference [27] also obtains an estimate for Csj (shown as 'B') based on the con- 
struction of a rational approximant which effectively coincides with the FTA ¥3 but 
using the physical value of = 92.4 MeV [20] instead of the value of fo in Eq. ([7]). 
Had they used fo, the result would have been lower, and would have agreed with the 
value we mentioned in the paragraph right after Eq. fllip . Therefore, our comments 
on the Fade Type found in that paragraph also apply to this determination in |27j . 
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Figure 1: Prediction for C^j in the large- A^^^ limit from the PA R^^f. [18], and the 

different PTAs discussed in the text and appearing in Table 1. For comparison we also 
show the estimate from Refs. [261 [271 [28] , which we label 'A', 'B' and 'C (resp.). 



Finally, one can get still another estimate for Cs7 from the PTA T2 in [27J by 
assuming that the ai mass in the large-A^^c limit is not approximated by the physical 
value in Eq. (ITTl) . but by a value which comes from the radiative pion decay saturated 
with the p and the ai. This value turns out to be ~ 998 MeV [28]. This lower 
number for the Oi mass is the reason for a higher value for Cgr than that obtained in 
[27], and is shown as 'C in Fig. 1. However, there is no compelling reason to associate 
this different mass of the ai with the large- A^c limit. In fact, our results show how 
similar values for Cs7 can be obtained with the physical masses of the mesons used 
for the poles. Morevover, one of the advantages of our method is that one can get 
a rough idea about the systematic error involved by looking at the dispersion of the 
values obtained. 

Of course, in the large- A'c limit Cgj does not run with scale whereas in the world 
at A'c = 3 it does. This is an additional source of systematic error in the result (fT3|) . 
However, phenomenological evidence as well as theoretical prejudice[23] suggests that 
a reasonable guess for this systematic error may be obtained by varying the scale in 
C%-j{ii) between the range 0.5 GeV < ^ 1 GeV (compare with Lio in ([7])). Using the 
running obtained in Ref. [25], this error turns out to be ~ 30 per cent, right in the 
ballpark expected for a typical l/N^. effect. This systematic error should be added to 
our large- A'c result in Eq. f[T^ in order to obtain an estimate for C^j{p ~ 0.7) in the 
real world. In this case, all the different results in Fig. 1 can be encompassed by this 
error. 

We would like to finish by recalling that PAs and PTAs are, in a way, two extreme 
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versions of a rational approximant. While in the latter all poles are fixed at the physical 
masses, in the former the poles are left free, and they are obtained by demanding that 
the expansion around = reproduces that of the original function to the highest 
possible order. Besides these two rational approximants, there are also the so-called 
Partial Pade Approximants [T8l [22] which, from a certain point of view, lie half way 
between PAs and PTAs. These Partial Pades are rational functions whose polynomial 
in the denominator has only some of the poles preassigned but the others are left free, 
to be determined by the usual matching conditions at = 0. Therefore, there is no 
reason why, in general, the poles of a Partial Pade should come out to be purely real§, 
unlike those of a PTA, which are of course real by construction. We have constructed 
seven of these Partial Pades, with a polynomial in the denominator up to fifth order 
in Q^. In some of the cases the poles were actually complex, as it was also the case 
of the PA P2 IIS]- However, the results obtained for Cgy are almost identical to those 
in Fig. 1, although with errors which are somewhat larger. This feature reinforces the 
stability of the result shown in Fig. 1, and gives us reassurance about the reliability of 
our result. Finally, we would like to mention that predicting 0{p^) LECs may also be 
another straightforward application of this method. 
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